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We study the effects of quasiparticle interactions in a quasi-two dimensional (quasi-2D), zero- 
temperature Bose-Einstein condensate of dipolar atoms, which can exhibit a roton-maxon feature in 
its quasiparticle spectrum. Our focus is the Beliaev damping process, in which a quasiparticle collides 
with the condensate and resonantly decays into a pair of quasiparticles. Remarkably, the rate for 
this process exhibits a highly non-trivial dependence on the quasiparticle momentum and the dipolar 
interaction strength. For weak interactions, the low energy phonons experience no damping, and the 
higher energy quasiparticles undergo anomalously weak damping. In contrast, the Beliaev damping 
rates become anomalously large for stronger dipolar interactions, as rotons become energetically 
accessible as final states. Further, we find a qualitative anisotropy in the damping rates when the 
dipoles are tilted off the axis of symmetry. Our study reveals the unconventional nature of Beliaev 
damping in dipolar condensates, and has important implications for ongoing studies of equilibrium 
and non-equilibrium dynamics in these systems. 


The quasiparticle picture of fluctuations and excited 
states in condensed matter systems is a fundamental 
modern paradigm. Early investigations in this direc¬ 
tion focused on superfluid ^He, which hosts very low en¬ 
ergy quasiparticles at intermediate wave vectors, termed 
“rotons” m- Rotons were first observed in neutron 
scattering experiments with “^He [n sHg, and are now 
understood to emerge in strongly interacting superfluids 
due to strong, longer-range two-body correlations [7HS]. 
Bose-Einstein condensates (BECs) of atoms with large 
magnetic dipole moments, such as Cr, Er, or Dy, are 
unique in that they are predicted to support roton quasi¬ 
particles when confined to highly oblate, quasi-two di¬ 
mensional (quasi-2D) geometries, despite remaining ex¬ 
tremely dilute and weakly interacting compared to super¬ 
fluid ^He [TUHTB] . Thus, mean-held theories typically pro¬ 
vide good descriptions of these systems [Ml [15] , despite 
their treatment of quasiparticles as free, non-interacting 
excitations. Here, we systematically step beyond the 
mean-held approximation, and study the effect of quasi¬ 
particle interactions on the damping of collective excita¬ 
tions in quasi-2D dipolar condensates, hnding non-trivial 
effects beyond the free quasiparticle picture. 

In 1958, Beliaev hrst presented a theory of the Bose- 
condensed state that includes quasiparticle interactions, 
showing how they manifest as effective condensate- 
mediated processes [laiiTi. An important consequence 
of such interactions is the damping of quasiparticle mo¬ 
tion, resulting in hnite lifetimes for collective condensate 
excitations. Beliaev specialized to the case of isotropic, 
short-range (contact) interactions, which is relevant for 
alkali atom condensates [TH]. A number of subsequent 
works have following along these lines [T^l - El] . and there 
is notable agreement with experimental work |251 I2hj . 
However, despite a growing interest in the experimen¬ 
tal study of quantum many-body physics with dipolar 
atoms P7H55] and polar molecules [51H55| . a systematic 


theoretical understanding of beyond mean-held effects, 
such as quasiparticle damping, is lacking for these sys¬ 
tems. 

In this Letter, we present a theory describing these ef¬ 
fects in a quasi-2D dipolar BEC, and hnd a number of 
striking results. When the dipolar interactions are weak, 
the damping rates are anomalously small, being signih- 
cantly less than those of a gas with purely contact inter¬ 
actions of equal strength. In contrast, when the dipolar 
interactions are stronger and rotons begin to emerge in 
the quasiparticle spectrum, the Beliaev damping rates ac¬ 
quire anomalously large values, though the rotons them¬ 
selves remain undamped. For all interaction strengths, 
the low energy phonon modes are immune to damping 
[39j . Additionally, the dipolar interactions can be made 
strongly anisotropic in the quasi-2D geometry [40] . In 
this case, the Beliaev damping rates acquire qualitatively 
different character depending on the direction of quasi¬ 
particle propagation; this feature has no analog in con¬ 
ventional superfluids. Our results mark an important 
step towards understanding the physics of dipolar con¬ 
densates beyond the mean-field approximation, and have 
important implications for both the equilibrium and non¬ 
equilibrium properties of these novel superfluids. 

In the grand canonical ensemble, the Bose gas Hamil¬ 
tonian is 

H = J drif^r) + U{r)- tp{r) 

-f - y dr y dr'?/;^(r)'!^^(r')H(r — r')^/)(r')7/)(r). (I) 

Here, m is the atomic mass, U{r) is the external po¬ 
tential, pL is the chemical potential of the gas, and '0(r) 
('0’^(r)) is the Bose annihilation (creation) operator. For 
fully polarized dipoles with dipole moments d, the two- 
body interaction potential is H(r) = d^(l — 3cos^ 9)/\r\^, 
where 9 is the angle between r and d. 
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At ultracold temperatures, the dilute Bose gas can be 
described by a mean-field theory with a condensate or¬ 
der parameter 4){y) = (^/>(r)), which evolves under the 
equation of motion, 

ih^{v) = +Uir)- f?j 0(r) 

-I- / dr'T^(r — r')('5^(r')d'(r')d'(r)). (2) 



Under the decomposition '&(r) = where cp{r) 

annihilates non-condensed atoms, (r')^(r)) ~ 

n(r')(^(r)-|-h(r',r)(/)(r'), where n(r) = |0(r)p-|-h(r, r) is 
the total density of the gas and h(r',r) = 
is the non-condensate density matrix. We work in the 
Popov approximation, and omit the anomalous density 
matrix m(r', r) = {(p{r')ip{r)) from the theory [41]. In the 
perturbative framework we employ, the Beliaev damping 
rates are insensitive to this approximation EOlllI]. 

Small amplitude condensate oscillations can be mod¬ 
eled as perturbations S(j){r) about the stationary state 
of Eq. denoted 0o(r). We obtain equations of mo¬ 
tion for these condensate oscillations by inserting (j){r) = 
(/)o(r) + S(l){r) into Eq. ([^ and linearizing about 5. If the 
couplings between S(j){r) and the non-condensate density 
n are ignored, this procedure reproduces the Bogoliubov 
free-quasiparticle description of small amplitude conden¬ 
sate oscillations. This description, however, is inadequate 
to describe the damping of condensate oscillations. To 
correct this, we couple the condensate oscillations, which 
take the form of Bogoliubov quasiparticles, to the non¬ 
condensate atoms perturbatively in 6, following the pro¬ 
cedures of Refs. [inumiMiiia- We obtain eigenfrequen- 
cies oj' = uj + Slu, where uj are the bare (non-interacting) 
quasiparticle frequencies and 6uj are frequency shifts that 
arise due to quasiparticle interactions. The imaginary 
part of Su! corresponds to a damping rate for conden¬ 
sate oscillations. At T = 0, this is a Beliaev process, 
which involves the resonant decay of a quasiparticle into a 
pair of quasiparticles under the constraints of energy and 
momentum conservation [I6lll7|. The relevant damping 
rate is thus 7 b = Im[(5a;]7’=o- This perturbative scheme 
remains valid for large damping rates, as long as the non¬ 
condensate density remains small. 

We restrict our study to the quasi-2D regime, where 
the atoms are free to move in-plane but are tightly 
confined in the axial direction by an external potential 
C7(r) = mw^z^/2. If is the dominant energy scale 
in the system, to a good approximation all atoms oc¬ 
cupy the single-particle ground state in the ^-direction, 
xiz) = where U = ^hjraiOz- An 

effective quasi-2D theory is obtained by separating all 
bosonic fields into this axial wave function and inte¬ 
grating the 2 -coordinate from the theory [13]. Below, 
we rescale all lengths in units of Iz and all energies in 
units of hu)z- The condensate order parameter becomes 


FIG. 1: (color online). Energetics of Beliaev damping for a 
quasiparticle with momentum (a) p4 = 0.95;, a maxon (black 
-f), and (b) p4 = 2.25;, in the free-particle part of the spec¬ 
trum (black circle), in a quasi-2D dipolar condensate with 
nogd = 1.7 and a = 0. The vertical panels show the quasi¬ 
particle spectrum (red lines). The dark blue lines show the 
manifold of decay channels allowed by energy and momentum 
conservation. The light blue lines show the allowed momenta 
only. 


4>{y) = y^noxiz) where ng is the uniform areal conden¬ 
sate density, and the condensate oscillations take the 
form 5<()(r) = x(^) 

where p and p are in-plane spatial and momentum co¬ 
ordinates, respectively. The coefficients Up and Up are 
the Bogoliubov quasiparticle amplitudes, given by Up = 
^£p/2bjp -b 1 and Vp = -sgn[U(p)]y'ep/2a;p - 1, where 
ffp = p^ 12 -b no 5 dU(p). The bare quasiparticle spectrum 
is 


(3) 

where gd = y/^d'^/3 is the quasi-2D dipolar interac¬ 
tion strength and V (p) = Fj_ (p) cos^ a -b Ey (p) sin^ a 
is the quasi-2D momentum-space dipolar interaction po¬ 
tential, with Ej_(p) = 2 — S^pe^ erfc[p] and E||(p) = 
— l+3x^{py/p)eP^eTic[p]. Here, erfc[p] is the complimen¬ 
tary error function and a is the polarization tilt angle 
between d = d{z cos a + y sin a) and the 2 -axis. The Be¬ 
liaev damping rate for a quasiparticle with momentum p 
is found to be 

9 _ 

= X E (4) 

kq 


where = A^q -b Aqi^., and A^q has matrix elements 

“^kq = Up (v (k) «nq -b 4uq) -b U(k -b qXwq 

+ (4^(k)(wk< + <<) + l^(k + q)uki^q) (5p.k+q- 

(5) 


We take the thermodynamic limit, and evaluate Eq. Q 
numerically. 

We first consider a quasi-2D dipolar condensate that 
is polarized perpendicular to the 2D plane (a = 0). In 
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this case, an expansion of the small-momentum, phonon 
part of the quasiparticle spectrum gives Wp ~ Cdp(l — 
/?>2p+ ...), where Cd = is the phonon speed. 

This downward curvature prohibits the Beliaev damping 
of phonons, due to the impossibility of simultaneous en¬ 
ergy and momentum conservation. Thus, phonons do 
not damp in quasi-2D dipolar condensates |39j . This 
is in contrast to quasi-2D condensates with repulsive, 
isotropic contact interactions, which host quasiparticle 
spectra with upward curvature at small momenta, re¬ 
sulting in Beliaev damping rates oc at small p |44j . 

At larger momenta, a local “roton” minimum with an 
energy gap Ai- develops in the quasiparticle spectrum for 
dipolar interaction strengths noPd ^ 1.15, which ulti¬ 
mately softens to Ar = 0 at a momentum p,. ~ 1.62 when 
noPd — 1.72. This is accompanied by a local “maxon” 
maximum at p ~ 0.74. An example roton-maxon spec¬ 
trum for riopd = 1.7 is shown in the vertical panels of 
Fig. and by the red curve in Fig. [^d). 

As the roton minimum develops, the density of quasi¬ 
particle states grows significantly. Near the minimum, 
the spectrum can be expanded about p ^ Pi- to give 
Wp ~ Ar -I- (p — pr)^/2mr where irir is the effective ro¬ 
ton mass. The density of states near the roton minimum 
is thus Pr{uj) = 27rmr(l + Pr/\/2mr(w — Wr)). The di¬ 
vergence of this expression at oj = u>r contributes to an 
anomalously large density of states in this vicinity. It 
is instructive to note that the expression for the Beliaev 
damping rate in Eq. Q is reminiscent of Fermi’s Golden 
Rule, which describes the scattering of a quantum state 
into other final states at a rate proportional to the den¬ 
sity of available final states. Indeed, the evaluation of 
the Dirac-delta function in Eq. 0 produces a factor re¬ 
sembling the density of final quasiparticle states; we thus 
expect large damping rates for quasiparticles that can 
decay into rotons. 

In Fig. we illustrate the manifold of available fi¬ 
nal quasiparticle states for no5d = I-!") which supports 
a prominent roton-maxon feature. In panel (a), we con¬ 
sider a quasiparticle with momentum p = 0.9x (shown 
by the black -|- sign in the Px-Py plane), which is in the 
maxon part of the spectrum. Whenever the maxon en¬ 
ergy exceeds 2Ar, it is energetically possible to decay 
into a pair of rotons. The blue lines, which show the 
energy and momenta of the available final quasiparti¬ 
cle states, are centered about the roton minima in the 
+y and —y directions. This indicates that maxons un¬ 
dergo Beliaev damping by decaying into a pair of nearly 
counter-propagating rotons that travel transverse to the 
initial quasiparticle direction. In panel (b) of Fig. we 
consider a quasiparticle with momentum p = 2.2i: (black 
circle), which is in the higher energy, free particle-like 
part of the spectrum. A number of final states are avail¬ 
able to these quasiparticles; they can decay into rotons, 
maxons, and phonons. In the latter process, the quasi¬ 
particle “sheds” low-energy phonons and loses a corre¬ 





FIG. 2: (color online), (a) Beliaev damping rate for weak 
dipolar interactions, where no roton is present in the quasi¬ 
particle spectrum; the corresponding spectra are shown in 
(c). (b) Rate for stronger dipolar interactions, where a roton 
is present; the corresponding spectra are shown in (d). 


spondingly small amount of energy and momentum. In 
the former processes, many combinations are final states 
are possible. Roton-maxon pairs can be produced, or 
pairs of forward-propagating quasiparticles; the momenta 
of these final states are shown by the detached blue loop 
in Fig. [^b). 

We plot the Beliaev damping rates for a range of dipo¬ 
lar interaction strengths in Fig. the rates are scaled 
by the axial trap frequency Panel (a) shows rates 
for quasiparticle spectra that lack roton-maxon features 
(shown in panel (c)). The downward curvature of the 
quasiparticle spectrum forbids phonon damping below 
a critical momentum Pcrit; this is apparent in all cases 
shown. For small p, pcrit = and Pcrit ~ Cd 

for larger p [31]. Additionally, we see evidence that as 
nopd increases, the downward curvature of the quasi¬ 
particle spectrum becomes more pronounced, and Pcrit 
increases correspondingly. As p —)■ Pcrit from above, 
the damping rate becomes anomalously large. The only 
available mechanism for Beliaev damping in this small 
range of momenta near pcrit is the shedding of low energy 
phonons. We attribute these anomalously large damping 
rates to the unique curvature of the quasiparticle spec¬ 
trum, which produces a factor resembling a large density 
of phonon states in the evaluation of Eq Q. 

Note that for UQPd = 0.1, the Beliaev damping rates 
are very small, remaining much less than ujz in the range 
of p shown. These rates are nearly an order of magnitude 
smaller than those of a quasi-2D non-dipolar BEG with 
an equivalent chemical potential. As n^gd increases, the 
damping rates increase significantly across the range of 
p, which we expect due to the proportionality 7 b oc g^. 
The rates at larger p become more comparable to those 
of a system with purely contact interactions as n^gd —>■ 1- 

Panel (b) of Fig. shows Beliaev damping rates for 
larger dipolar interaction strengths, which support spec- 
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tra with pronounced roton-maxon features (shown in 
panel (d)). We consider two distinct cases; for n^gd = 
1.5, the maxon energy is less than 2Ar (blue curve), and 
for riQgd = 1.7, the maxon energy is greater than 2A,. 
(red curve). In the former case, it is energetically forbid¬ 
den for a maxon to decay into a pair of rotons. Thus, all 
low-energy quasiparticles (phonons, maxons, and rotons) 
remain undamped, and Beliaev damping only occurs for 
p > 2. For the latter case, a maxon can damp into a pair 
of transverse, counter-propagating rotons. Notice that 
the red curve in panel (b) is separated into three distinct 
parts. The two left-most parts correspond to Beliaev 
damping into roton pairs only. These damping rates are 
anomalously large, achieving values well over lOOa;^ for 
some values of p. This is due to the large density of 
states near the roton minimum. Additionally, the Beli¬ 
aev damping rate vanishes for a range of p near the roton 
minimum, reflecting the fact that rotons are undamped, 
due to their anomalously low energy and large momenta. 
The black -|- sign and circle show the Beliaev damping 
rates for quasiparticles with p = 0.9 and p = 2.2 respec¬ 
tively, corresponding to the discussion of Fig. Though 
the non-condensate density grows as A^ softens [U, it 
remains dilute for the cases we consider here. We thus 
expect our perturbation theory to remain valid for these 
large damping rates. 

By tilting the external polarizing field off axis (a 7 ^ 0), 
the dipolar interactions can be made strongly anisotropic. 
It is predicted that anisotropic dipolar interactions will 
produce a quasiparticle spectrum with correspondingly 
strong anisotropies, supporting rotons for only a narrow 
range of propagation directions m- Such anisotropic 
spectra have important consequences for the damping of 
quasiparticles in these systems. 

We plot the quasiparticle spectra and Beliaev damp¬ 
ing rates for a condensate with n^ga = 1-3 and a tilt 
angle a = tt/S in Figs, [^b) and l^c), respectively. In 
this case, the spectrum for quasiparticles propagating in 
the a:-direction (T, red line) exhibits roton-maxon char¬ 
acter, while the spectrum in the y-direction (||, blue line) 
does not. Above, we noted that maxons can only un¬ 
dergo Beliaev damping by decaying into a pair of nearly 
counter propagating rotons when a = 0. Here, no ro¬ 
tons exist in the transverse direction, and maxons are 
consequently undamped despite the fact that the maxon 
energy exceeds 2Ar. Quasiparticles propagating in the 
cc-direction begin to damp near p = 1.7i, shown by the 
black circle(s) in Fig. The momenta of the available 
final states are shown by the red line in panel (a), and 
the corresponding damping rates are shown in panel (c). 
The onset of damping is due to the shedding of phonons 
near this momentum. Interestingly, the damping rate is 
not anomalously large near this onset, unlike the a = 0 
case; this is due to the anisotropy of the spectrum, which 
skews its curvature unfavorably. 

Although no roton-maxon feature exists in the y- 



p 



FIG. 3: (color online). Beliaev damping of a quasi-2D dipo¬ 
lar condensate with nopd = 1.3, a = tt/S, and 9 — 0 (dipoles 
tilted in the t/-direction). (a) Manifold of allowed decay chan¬ 
nels for a quasiparticle with momentum p = 0.7y (black -I-) 
shown by the blue lines, and momentum p = 1.7x (black cir¬ 
cle) shown by the red lines, (b) Spectrum for quasiparticles 
propagating in the j/-direction (blue line) and the ^-direction 
(red line), (c) The corresponding Beliaev damping rates. 


direction, quasiparticles propagating in this direction can 
damp by decaying into transverse roton pairs, as illus¬ 
trated by the blue lines in Fig. [^a), which show the 
momenta of the available final states for a quasiparticle 
with p = \.7x (shown by the black -|- sign). The damp¬ 
ing rates for this process are shown by the two left-most 
blue line segments in panel (c). For larger momenta, the 
quasiparticles begin to shed phonons in the j/-direction. 
The damping rates for these momenta are shown by the 
right-most blue line in panel (c). Thus, the anisotropic 
dipolar interactions not only lead to anisotropic damping 
rates, but rather qualitatively different damping mecha¬ 
nisms depending on the direction of quasiparticle propa¬ 
gation. 

Our predictions have important consequences for on¬ 
going experiments with ultracold dipolar atoms. For 
example, in experiments measuring the dynamic struc¬ 
ture factor S'(p,a;) of the condensate via, for exam¬ 
ple, optical Bragg scattering [i^, these rates will ap¬ 
pear as spectral widths |3S]. Further, our results can 
be used to predict the short-time non-equilibrium dy¬ 
namics of these systems, as the Beliaev mechanism is 
responsible for the redistribution of quasiparticles near 
T = 0. Take, for example, the anisotropic case dis¬ 
cussed above. If an oblate dipolar condensate is pre¬ 
pared with nogd = 1.3 and a = tt/S, and modes with 
p = 0.5a; are excited, they should undergo coherent dy¬ 
namics for long times. On the other hand, the excitation 
of modes with p = 0.5y will result in the nearly imme¬ 
diate redistribution of energy into transverse rotons. In 
this sense, the anisotropic Beliaev damping should re- 
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suit in strongly anisotropic relaxation dynamics. Exper¬ 
imentally, the limit of a deep roton {n^gd = 1-7) can 
be achieved, for example, with ^^'^Dy [29] in an oblate 
trap with axial frequency = 27r x 10^ Hz and a mean 
3D density han ^ 3 x lO^'^cm”^. For strongly dipolar 
molecules, much smaller densities are required. 
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